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Abstract 

It is suggested and demonstrated that two specific 2-dimensional 
correlation patterns, fixed-to-arbitrary bin and neighboring bin corre- 
lation patterns, are efficient for identifying various random multiplica- 
tive cascade processes. A possible application of these two correlation 
patterns to single event analysis in Relativistic Heavy Ion Collider 
experiments is discussed. 
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It is well-known that nonlinear phenomena are very popular in nature. Once a set of anomalous 
scaling law of probability- or "mass" -moments with respect to spatial size [|l| is observed in experi- 
ment, we recognize that there is self-similar multifractal structure. However, it is difficult to judge 
what kind of mechanism most likely causes the results since the conventional measures have not 
provided sufficient information on underlying dynamics. This hinders us from further understanding 
and studying nonlinear physics. A useful measure for this purpose should not contain too much 
detail, but enough to capture the essence of production mechanism. In this letter, as a first try 
along this direction, we will suggest two specific correlation patterns which manifest distinguishable 
characters of various random multiplicative cascade processes. They will be helpful in exploring the 
origin of nonlinear and correlation related phenomena. 

A random multiplicative cascade process is the simplest example of having a well-defined mul- 
tifractal structure and has been used extensively in various fields. For example, energy dissipation 
in fully developed turbulence |^ and successive branchings in QCD parton-shower ^ are both 
this kind of process and accordingly have multifractal structure. The main idea of such a process 
is simply a series of self-similar random cascades in spatial partitions. It is generally described as 
follows: at the first step of the cascade, a given initial interval with length A and unit probability 
Po = 1 is split into r subdivisions ("bins"). Usually, r is an integer with r > 2. The probability for 
a bin to be split is determined by a weight w which is distributed according to a splitting function 
p{wi, W2, ■ ■ ■ , Wr). Then each bin obtained from the first step of the cascade is again independently 
split into r bins and so on. After u = 1,2, J steps, the final number of bins is t'^, and the probabil- 
ity in a specific bin is z/ products of w over all previous generations: pj^^ = Wij^W2j2 ■ ■ ■ Wyj^, where 
= 1, 2, r'' (/i = 1, 2, z/). The distribution of pj^ obtained in this way fiuctuates violently from 
bin to bin, cf. Fig. 6 in [^. 

By varying the splitting number and/or changing the splitting function, all kinds of random 
multiplicative cascade processes can be constructed. Especially, an additional physical constraint 
on the splitting function, like probability conservation at each step of splitting, will cause stronger 
correlations among the r split bins than those without it. The smaller the splitting number r is, 

2 



the stronger are the correlations among the r bins. It is therefore necessary to investigate n-point 
correlations in order to probe the differences among the various random cascade processes. Let us 
start from the simplest 2-point ones. 

Why are conventional measurements for 2-point-correlation moments @] insufficient? This is 
because they are usually defined at fixed correlation length with both horizontal average over all 
bins in an event and longitudinal average over all events in a sample. The fixed correlation length 
makes a comparison between correlations with different scales impossible, and the horizontal average 
smooths out the differences between different pairs of bins. For these measurements, various random 
multiplicative processes have the same kind of power law behavior with diminishing correlation 
length, similar to the scale invariance of "mass" moments. They have little help in the identification 
of production mechanism. Hence, 2-point correlation moments without horizontal averaging, but 
with various correlation lengths, are necessary to be presented simultaneously. It should be a pattern 
constructed by all kinds of 2-point correlation moments. 

A straightforward pattern of such a kind is the three dimensional pattern of 2-point correlation 
cumulants 0], 

CkuK2 = (Inm lnpi^2)c = {\npK, Inpj^^) - {\npK,){\npK2), (1) 

where Ki and K2 are the coordinates of two points. However, in this pattern the differences among 
various random multiplicative cascade processes are hidden in complicated background and hard to 
be observed, cf. Fig.l, where the patterns for three different models, the a model without probability 
conservation |^, the p model with probability conservation ||^ and the c model with both probability 
conservation and random probabihty partition at each step of the splitting 0, are presented. The 
pattern manifests the common characteristics of random multiplicative cascade processes, such as 
symmetry at each step of splitting, the largest correlations of any bin with itself and so on, but 
is insensitive to the basic model assumptions. Thus, it is better to present the correlations of 
some specific two bins in a two dimensional pattern so that the differences in various processes are 
manifested. 
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A simple way to do so is to fix one bin Ki and vary the left one K2. It is a fixed-to- arbitrary 
bin correlation pattern. It contains correlation information at various lengths from the fixed bin to 
all other bins, in which r — 1 bins are separated from the fixed one at the last cascade step and 

— T bins are separated successively from the fixed one before the last step of cascading. Another 
simple way is to fix the distance between the two bins and change the position of them. Since the 
correlations among r bins at each step of splitting are important information, the distance between 
two bins are better to be as close to each other as possible, i.e. Ki = K and K2 = K + 1, where 
K = 1,2, ...,r'^ — 1. This is a neighboring bin correlation pattern. These two correlation patterns 
provide useful information on cascade mechanism and filter out the common characteristics of random 
multiplicative cascade processes. 

In the following, as a support of above arguments, the correlation cumulants for three different 
random cascade processes, the a, the p and the c models, are first derived in order to show how the 
different appearances of these cumulants are caused by the model assumptions. Then it is argued 
and demonstrated why the measurements for these two patterns provide information on the model 
assumptions. Finally, the extension and application of the patterns are discussed. 

These three models are examples of the random binary cascade processes with symmetric splitting 
functions ^(^1,^2) = 'p{.'^2,Wi). In this case, Eq.(l) can be written simply as 



are respectively the so called "same-lineage" cumulant, which is the correlation of the common 
ancestor of Ki and K2, and the "splitting" cumulant, which measures the correlation between the 
two parts split first from their common ancestor; 



(2) 



where 62 = J — Z]f=i ' ' ' k^ is the ultrametric distance, which is a measure of how many 
generations one has to move up before a common ancestor is found; 



A = d^Q[\i,\2]/d\l\x=o 



B = d^Q[\l,\2]/d\idX2\x=0 



(3) 
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is the branching generating function (BGF). 

For the a model, the Wi or W2 at each step of sphtting is two possible numbers, 1 + a and 1 — a, 
with equal probability for each of them. Hence, its splitting function P] can be written as: 



p{wi,W2) 



5[ wi-{l+a) ]+5[ wi-{l-a) 



5[ W2-{l+a) ]+5[ W2-{l-a) 



(5) 



where a is a fixed parameter in the region [-1, 1]. Inserting this function into Eq.(4) and then Eq.(4) 
into Eq.(3), the "same-lineage" and the "splitting" cumulants can be easily derived: 



In 



l + a 
1 — a 



B^ = 0. 



(6) 



Ba = means that there is no correlation between the two parts for every splitting. This is due to 
the independence of Wi and W2 assumed in the model. 

For the p model, the summation of Wi and W2 keeps to be 2, with Wi = 1 ± [3 and W2 = 1 =F /3 
as the consequence of probability conservation, here /5 is a fixed parameter in the region [-1, 1]. Its 
splitting function H] is 



p{Wi,W2) 



5{w,-il+l3) +6{w,-il-l3) 



6{wi + W2~ 2), 



(7) 



and its "same-lineage" and "splitting" cumulants are 



In- 



1-/5 



In 



1 + 



1-/5 



Due to probability conservation, the "splitting" cumulant B has sign opposite to the "same-lineage" 
cumulant A. Moreover, since /? is a unique fixed number, the derivatives of BGF with respect to one 
of the two A's twice are equivalent to those with respect to both of them. This is why the absolute 
value of the "same-lineage" cumulant A is equal to that of the "splitting" cumulant B. 

For the third c model, Wi + W2 = l similar to the p model. Distinct from the p and the a models, 
the weight w is allowed to be a random number e.g. Wi = {1 + 7r)/2 and W2 = {1 — 7^)/2, here r is 
a random number uniformly distributed in the interval [—1, 1] and 7 is a fixed model parameter in 
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the region [0, 1]. Obviously, this model is more flexible and closer to the real system with probability 
conservation. The splitting function is: 



p{wi,W2) = - 

7L V 2 

and its "same-lineage" and "splitting" cumulants are: 



— — 0\wi — 



5{wx^W2-l). (9) 



1 1 r 

Ac^- I ' {Inwifdwi I ^ Inwidwi 

7 7^ Li V 

1 1 r /-^ 

Be = — Inwi ln(l — wi)dwi ^ Inwidwi 
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The "sphtting" cumulant is nonzero due to probability conservation. The absolute values of the 
"same- lineage" and the "splitting" cumulants are unequal since Wi and W2 are no longer fixed to a 
particular number, unlike the p model. 

Prom the above analytic derivation, we see that the three random cascade models are very different 
in the absolute values of their "splitting" cumulants and in the relative values of these cumulants 
to the "same lineage" cumulants. These diff'erences arc due to the basic model assumption and 
independent of the model parameters. Therefore, a useful 2-point correlation pattern should contain 
information on the "splitting" cumulant and its relation to the "same hneage" cumulant. This 
requirement turns out to be two possible combinations of the coefficients A and B in Eq.(2): (1) 
Cki,K2 ^ B, i.e. J - c?2 = and 1 - ^^2,0 = 1; (2) Cki,K2 ^A + B, i.e. J - d2 = I and 1 - 5^2,0 = 1; 
for the first combination, d2 — J and ^2 7^ means that common ancestor has to be at the initial 
interval if J > 0. The fixed-to-arbitrary bin correlation pattern happens to have 2^~^ • (2 — 1) pairs of 
this kind of bins. The neighboring bin correlation pattern also has a pair of this kind of bins at center 
K — 2"^"^. For the second combination, ci2 = — 1 and d2 ^ implies that common ancestor locates 
at the first cascade step. In the fixed-to-arbitrary bin correlation pattern, 2'^~^ • (2 — 1) pairs of bins 
have common ancestor at the first cascade step. The neighboring bin correlation pattern again has 
two pairs of such a kind of bins at K — 2*^"^ and K — 2'^"^ • 3 respectively. Hence the suggested 
two correlation patterns are useful measures in identifing the above mentioned three random cascade 
processes. 



We first plot the fixed-to-arbitrary bin correlation pattern for the a, the p and the c models 
in Figs. 2(a), (b) and (c), respectively, with the model parameters a — 0.4, /3 = 0.4, 7 = 0.8, 
J = 6. The ordinate is the correlation strength Ci^k2'i the abscissa is the bin K2 which varies as 
K2 — 2, 3, 64, and the fixed bin is located at Ki — 1. Since the larger the distance between the two 
bins is, the earlier the two bins are separated in cascade process, the correlation strength decreases 
as the arbitrary bin K2 goes further away from the fixed one Ki. The number of points with the 
same height or same correlation strength increases as 2"^^ ■ (2 — 1) with v = 1, 2, J. The number 
2 comes from the binary cascade. Generally, it is t"'^ • (t — 1). So by counting the number of points 
in each correlation strength, the splitting number of the cascade process can be estimated. 

The possible correlation strengths in these three models are different: (1) the weakest one in the 
a model is zero; (2) there are both zero and negative ones in the p model; (3) negative ones, but no 
zero ones, are allowed in the c model. These differences come from the different model assumptions. 
The absence of negative correlations imphes independent splitting in the model, B — 0, such as is the 
case for the a model. If there are both negative and zero correlations in the pattern, the positive one 
from the "same-lineage" cumulant A must be equal to the negative one from the "splitting" cumulant 
B. So that S < for the weakest one and A + B — ior the following one, as for the p model. 
Having negative, but no zero, correlation, means that the "same-lineage" cumulant is different from 
the "splitting" cumulant, or A + B ^ and B < 0, as for the c model. 

The neighboring bin correlation pattern Ck,k+i with K — 1,2, 2"^ — 1 can show us the same 
thing. In Fig. 3, the patterns for the three models are presented. Here, the strongest correlations are 
of those two neighboring bins split at the last step of cascade, i.e., K = 1,3, ...,2'' — 2°; the next 
ones happen to be of those split before the last step of cascade, i.e., K — 2,6, 10, 2-^ — 2^; and so 
on. The number of points at a certain correlation strength is determined by 2^"^, or t'^'~^ in general, 
with u = J, J — 1, ...,1. There are only positive and zero correlation spectra in the a model. In the p 
model, the correlation spectrum ranges over positive, zero and negative values, while for the c model, 
only positive and negative correlations are allowed with no zero ones. 

Therefore, from the measurements of these two correlation patterns, we can find out by what 
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kind of random multiplicative cascade process the observed system are produced. This is a first step 
toward the revelation of underlying mechanism of multifractal phenomena. The extension of these 
two patterns to the mixture of different random cascade processes and other mechanisms, which 
can lead to multifractal structure, is obviously the next step. Their corresponding patterns will 
not be so regular as those of the above mentioned three processes. For the further identification 
of the production mechanism in general, other measures may be needed in addition to these two. 
Nevertheless, having established these two patterns for some known mechanisms, it is interesting 
to compare the patterns measured in practice with them and judge if there are similar generating 
mechanism. 

A possible application of these two correlation patterns is in the single event analysis of Relativistic 
Heavy Ion Collider (RHIC) experiments p. In these experiments, the probability of final state 
particles falling into a particular bin in rapidity or pseudo-rapidity can be well estimated by 0, 
~ ukJN, where is the number of final state particles falling into the Ki-th bin and is 
the total number of particles in the event, which is typically several thousand [|T^. At the current 
experimental resolution in rapidity, there is almost no empty bin. In this case, it is much better to 
investigate directly the pattern of a single event without averaging over all events in the sample. The 
2-point correlation cumulant for a single event can be defined as: 

CS,K2 = l^P^i 1^P^2 - (lnpi^i)(lnp/^2), (12) 

its average over all events in a sample equals to the 2-point correlation cumulant Cki.K2 of Eq.(l). 
The whole event sample can then be classified by the measured correlation patterns of single event. 
Some events which have undergone a phase transition to Quark Gluon Plasma (QGP) will most 
likely have very special patterns, since during this transition, all kinds of correlations, the long as 
well as the short ones, will change dramatically. 

In this letter, we investigate two specific 2-dimensional correlation patterns: the fixed-to-arbitrary 
bin and the neighboring bin correlation patterns. The former shows correlations at various lengths, 
and the latter has short-range correlations at different positions. They are sensitive to the underlying 
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mechanisms and work equally well in identifying various random multiplicative cascade processes. 
As an example, their application to the analysis of single event of current Relativistic Heavy Ion 
Collider (RHIC) experiments has been discussed. For more-than-2-point correlations, we can also find 
some useful correlation patterns in a similar way. The measurements for all of them will hopefully 
open a new world in the research of nonlinear and correlation related physics. 

We are grateful for the helps of Dr. Jimmy MacNaughton in the presentation of the paper and 
the good suggestions of Dr. Nu Xu. This work is supported in part by the SFC under project 
90103019; Cross-century Talent Foundation of National Education Committee of China under Grant 
No.l [1998]. 
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Figure Captions 



Fig.l The correlation cumulants Cki,K2 oi (a) the a model, (b) the p model and the 
(c) c model. 

Fig. 2 The fixcd-to-arbitrary bin correlation patterns Ci^Xi of (a) the a model, (b) the 
p model and (c) the c model, where the ordinate is correlation strength Ci^k2 
and the abscissa is K2 — 2, 3, • • • , 64. The sohd hnes in the figures indicate 
zero correlation strength. 

Fig. 3 The neighboring bin correlation patterns Ck,k+i of (a) the a model, (6) the p 
model and (c) the c model, where the ordinate is correlation strength Ck,k+i 
and the abscissa is K — 1, 2, • • • , 63. The solid lines in the figures indicate zero 
correlation strength. 
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